Solitons in one-dimensional nonlinear Schrodinger lattices with a local inhomogeneity 
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In this paper we analyze the existence, stability, dynamical formation and mobility properties 
of localized solutions in a one-dimensional system described by the discrete nonlinear Schrodinger 
equation with a linear point defect. We consider both attractive and repulsive defects in a focusing 
lattice. Among our main findings are: a) the destabilization of the on-site mode centered at 
the defect in the repulsive case; b) the disappearance of localized modes in the vicinity of the 
defect due to saddle-node bifurcations for sufficiently strong defects of either type; c) the decrease 
of the amplitude formation threshold for attractive and its increase for repulsive defects; and d) 
the detailed elucidation as a function of initial speed and defect strength of the different regimes 
(trapping, trapping and reflection, pure reflection and pure transmission) of interaction of a moving 
localized mode with the defect. 
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I. INTRODUCTION 



The past few years have witnessed an explosion of in- 
terest in discrete models that has been summarized in 
a number of recent reviews This growth has been, 
to a large extent, motivated by numerous applications 
of nonlinear dynamical lattice models in areas as broad 
and diverse as the nonlinear optics of waveguide arrays 
0, the dynamics of Bose-Einstein condensates in peri- 
odic potentials [|[ , micro- mechanical models of cantilever 
arrays Q , or even simple models of the complex dynam- 
ics of the DNA double strand Arguably, the most 
prototypical model among the ones that emerge in these 
settings is the, so-called, discrete nonlinear Schrodinger 
equation (DNLS) @, 0. DNLS may arise as a direct 
model, as a tight binding approximation, or even as an 
envelope wave expansion: the DNLS is one of the most 
ubiquitous models in the nonlinear physics of dispersive, 
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discrete systems. 

Perhaps the first set of experimental investigations 
that generated an intense interest in DNLS type equa- 
tions was in the area of nonlinear optics and, in partic- 
ular, in fabricated AlGaAs waveguide arrays Q. In the 
latter setting a wide range of phenomena such as dis- 
crete diffraction, Peierls barriers (the energetic barrier 
that a wave needs to overcome to move over a lattice — 
sec details below), diffraction management (the periodic 
alternation of the diffraction coefficient) [jjllicj and gap 
solitons (structures localized due to nonlinearity in the 
gap of the underlying linear spectrum) (llj among others 
[12| were experimentally observed. These phenomena, in 
turn, led to a large increase also on the theoretical side of 
the number of studies addressing such effectively discrete 
media. 

A related area where DNLS, although it is not the 
prototypical model, it still yields useful predictions both 
about the existence and about the stability of nonlinear 
localized modes is that of optically induced lattices in 
photorefractive media such as Strontium Barium Nio- 
bate (SBN). Since the theoretical inception of such a 
possibility in Ref. |13l|. and its experimental realization 
in Refs. [Til [Tol Il6|. there has been an ever-expanding 
growth in the area of nonlinear waves and solitons in 
such periodic, predominantly two-dimensional, lattices. 
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A wide array of structures has been predicted and ex- 
perimentally observed in lattices induced with a self- 
focusing nonlinearity, including, e.g., discrete dipole [13], 
quadrupole [Hl| . necklace and other multi-pulse pat- 
terns (such as e.g., soliton stripes [III), discrete vortices 
pll ] , and rotary solitons [HJ • Such structures have a def- 
inite potential to be used as carriers and conduits for 
data transmission and processing, in the setting of all- 
optical communication schemes. A recent review of this 
direction can be found in Ref. [23| (see also Ref. [HI). 

Finally, yet another independent and completely dif- 
ferent physical setting where such considerations and 
structures are relevant is that of soft-condensed mat- 
ter physics, where droplets of the most recently discov- 
ered state of matter, namely of Bose-Einstein conden- 
sates (BECs) , may be trapped in an (egg-carton) optical 
lattice (OL) potential produced by counter-propagating 
laser beams in one, two or even all three directions [25j |. 
The field of BEC has also experienced a huge growth over 
the past few years, including the prediction and manifes- 
tation of modulational instabilities (i.e., the instability 
of spatially uniform states towards spatially modulated 
ones) [26], the observation of gap solitons [27j |. Landau- 
Zcner tunneling (tunneling between different bands of the 
periodic potential) [28| and Bloch oscillations (for matter 
waves subject to combined periodic and linear potentials) 
[29| among many other salient features; reviews of the 
theoretical and experimental findings in this area have 
also recently appeared in Refs. [30l l3ll]. 

While DNLS combines two important features of many 
physical lattice systems, namely nonlinearity and period- 
icity, yet another element which is often physically rele- 
vant and rather ubiquitous is disorder. Localized impu- 
rities are well-known in a variety of settings to introduce 
not only interesting wave scattering phenomena [32| , but 
also to create the possibility for the excitation of impu- 
rity modes, which are spatially localized oscillatory states 
at the impurity sites [33j . Physical applications of such 
phenomena arise, e.g., in superconductors 13411 . in the 
dynamics of the electron-phonon interactions [35| , in the 
propagation of light in dielectric super-lattices with em- 
bedded defect layers [36| or in defect modes arising in 
photonic crystals (37| . 

In the context of the DNLS, there have been a number 
of interesting studies in connection to the interplay of 
the localized modes with impurities. Some of the initial 
works were either at a quasi-continuum limit (where a 
variational approximation could also be implemented to 
examine this interplay) [38[ or at a more discrete level 
but with an impurity in the coupling [33 ] (see also in the 
latter setting the more recent studies of a waveguide bend 
poLlill] and the boundary defect case of Ref. [42|]). More 
recently the experimental investigations of Refs. [43|, |44j 
motivated the examinations of linear [45) ] and nonlinear 
[451I46T] defects in a DNLS context. In the photorefractive 
context, further recent experimental work has illustrated 
blocking effects to a probe beam from either bright or 
dark soliton beams in defocusing waveguide arrays [47l |. 



Our aim in the present work is to systematically ex- 
amine the properties of the focusing DNLS equation in 
the presence of both an attractive and a repulsive linear 
impurity. Our first aim is to present the full bifurcation 
diagram of the localized modes in the presence of the 
impurity and how it is drastically modified in compari- 
son to the case of the homogeneous lattice. The relevant 
bifurcations are quantified whenever possible even ana- 
lytically, in good agreement with our full numerical com- 
putations. A second problem that is examined is that 
of the threshold for the formation of solitary waves and 
how it is systematically affected by the presence of im- 
purities both in the repulsive and in the attractive case. 
This is motivated by the recent examination of the rel- 
evant threshold in the homogeneous lattice [48| and its 
connection with experiments in focusing [49| (and even 
defocusing [501 ]) waveguide arrays. Finally, in the same 
spirit as that of Ref. [451 ] , but for attractive and repulsive 
impurities, we systematically investigate the interaction 
of an incoming solitary wave with the localized impurity, 
identifying the main observed regimes as being trapping, 
reflection with trapping, pure reflection and pure trans- 
mission. 

This paper is organized as follows. In section II, we 
introduce the model. In Section III, we analyze the exis- 
tence and stability of localized excitations in a system de- 
scribed by the DNLS with the linear impurity. In Section 
IV, we examine the (energy /initial amplitude) threshold 
for soliton formation. In Section V we present our results 
related to the interaction of moving localized excitations 
with the impurity and, finally, in Section VI, we summa- 
rize our findings and present our conclusions. 

II. THE MODEL 

We consider a discrete system with a defect described 
by the DNLS equation as 

iipn + l\ipn\ 2 ^n + C(ip n+ i + ip n -i) + a n i/j n = 0, (1) 

where ip n is the complex field at site n (n = 1 . . . N) ; 7 
is the anharmonicity parameter, C the coupling constant 
and parameters a n allow for the existence of local, linear 
inhomogeneities. In this paper, we consider a single point 
defect, thus a n — a5 n ^ no , that can be positive (attractive 
impurity) or negative (repulsive impurity). In general, 
the presence of an on-site defect would affect the near- 
est neighbor coupling, and Eq. fl} should be modified 
to take this effect into account, as in Ref. [HI]. O n the 
other hand, this inhomogeneity in the coupling can be 
avoided using different techniques, for example, in non- 
linear waveguide arrays, changing slightly the separation 
between defect waveguide and its nearest neighbors, as it 
has been done in Ref. Q. In this work, we will assume 
that the coupling parameter C is independent on the site 
and positive. 

Upon renormalization of parameters, we consider 7 
1 (focusing case). Note that the defocusing case (7 < 
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0) can be reduced, under the staggering transformation 
tp n — > (— l) n ip n , to the previous one with opposite sign 
of the impurity a. Also, under the transformation ip n — > 
i/i„e 2 Eq. (fTJ) can be written in the standard form 



ilpn + l\^n\ 2 1pn + C Alp n + U n l\) n = 0, 



(2) 



where Aip n — ip n +i + ip n -i ~ 2^„ is the discrete Lapla- 
cian. Throughout this work, we use the form given by 
Eq. Q). 

The DNLS (JTJ) conserves two dynamical invariants, the 
Hamiltonian 



a > (attractive impurity): 



and 



V2C 



= 2C\ 



a n n ~ N 



( 



(8) 



(9) 



with an in-phase pattern (see bottom-right panel in 

Fig.rrj. 

a < (repulsive impurity): 



H = -J2l\^\ 4 +C(r n ^n+l+rn^n-l)+a n m 2 , (3) 0n = (-l)"0O 



with canonical variables q Tl = ^„ and p n = iip^i an< ^ the 
(squared L 2 ) norm or optical power 



(4) 



III. STATIONARY SOLUTIONS 

In order to study solitons in the system described by 
Eq. ([T]) , we aim to look for stationary solutions with fre- 
quency u. Thus, substituting 



(5) 



The stationary analog of Eq. (fTJ then reads 

- U>4> n + C{<p n+ i + </> n -l) + <i>\ + Ot n <t> n = 0. (6) 

Some of the properties of solitons are related to the ex- 
istence (or not) and properties of linear localized modes. 
These modes arise when an inhomogeneity appears and 
can be obtained from the linearized form (around the 
trivial solution <j> n — 0, V n) of Eq. ([6]). In this case, 
and considering an inhomogeneity located at the first site 
of the chain and with periodic boundary conditions, the 
problem reduces to solving the eigenvalue problem 
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that is a particular case of the eigenvalue problem studied 
in Ref. [53]. There it was shown that, if a ^ 0, the 
solution corresponds to TV — 1 extended modes and an 
impurity localized mode. Also, if N becomes large, the 
frequencies of extended modes are densely distributed in 
the interval fl £ [— 2C, 2C] and the localized mode can 
be approximated by 



and 



V2C 



-2CI 



■0 



-N 



0: 



or 
4C2 



+ 1 



(10) 



(11) 



with a staggered pattern (see bottom-left panel in Fig.[l|. 
In both cases 4>o is an arbitrary constant. In Fig. [1] we 
depict the linear mode spectrum as a function of the inho- 
mogeneity parameter a (top panel) and examples of the 
profiles of the ensuing localized modes (bottom panels). 

In order to explore the existence and stability of the 
nonlinear stationary states described by Eq. ([5]) , we have 
used the well-known technique based on the concept of 
continuation from the anti-continuum (AC) limit using a 
Newton- Raphson fixed point algorithm [53[ . Also, a stan- 
dard linear stability analysis of these stationary states 
has been performed, using the ansatz 

0n = [0soi + e(a n exp(Ai) + b„ exp(A*£)] exp(wt), (12) 

and solving the ensuing eigenvalue problem. <p so \ is the 
solution of Eq. ((6]) with frequency w, A is the linearization 
eigenvalue and A* its complex conjugate. Due to symme- 
tries of the system, the eigenvalues appear in quartets (if 
A is an eigenvalue, so are A*, —A and —A*). Furthermore, 
the U(l) invariance of the equation (the so-called phase 
or gauge invariance) leads to the existence of a pair of 
zero eigenvalues. If the remaining eigenvalues are imagi- 
nary, the state is linearly stable and, on the contrary, the 
presence of a eigenvalue with a nonzero real part implies 
instability. 

In the homogeneous lattice case of a = 0, fundamental 
stationary modes are well known to exist and be centered 
either on a lattice site or between two adjacent lattice 
sites [f| . The site-centered solitary waves are always sta- 
ble, while the inter-site centered ones are always unstable 

i- 

In order to study the effects of the inhomogeneity on 
the existence and properties of localized modes, we have 
performed a continuation from the homogeneous lattice 
case of a = 0. We found that, if a increases, (a > 0, 
attractive impurity case), the amplitude of the stable on- 
site mode decreases, while if a decreases (a < 0, repulsive 
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FIG. 1: Linear modes: the top panel shows the dispersion 
relation as function of impurity parameter a (notice the im- 
purity mode outside of the interval [— 2C, 2C]). The linear 
modes are normalized l<^n| 2 = l)i the impurity is lo- 

cated at n = 0, and periodic boundary conditions are consid- 
ered. The bottom panels depict examples of the profiles of 
the impurity modes. Bottom- left panel: profile for a = — 1 
(repulsive impurity) , and bottom-right panel profile for a — 1 
(attractive impurity). In all cases N = 200 and (7=1. 



impurity case), in general, the stable on-site soliton local- 
ized at the impurity merges with the unstable inter-site 
centered one localized between impurity and its neigh- 
boring site (beyond some critical value of |a|) and the 
resulting state becomes unstable. Notice that, at heart, 
the latter effect is a pitchfork bifurcation as the on-site 
mode collides with both the inter-site mode centered to 
its right, as well as with the one centered to its left. 

In Fig. [5] we show a typical bifurcation scenario where, 
for a fixed value of the frequency u> and the coupling pa- 
rameter C , we depict the mode power P corresponding 
to different on-site and inter-site localized modes as a 
function of impurity parameter a. If we denote as no the 
site of the impurity, when a > increases, we found that 
the unstable intersite soliton localized at n = n + 0.5 
disappears in a saddle-node bifurcation with the stable 
site soliton localized at n = no + 1. Also, if we continue 
this stable mode, when a decreases, and for a given value 



FIG. 2: Bifurcation diagram of stable (solid line) and unstable 
(dashed line) nonlinear modes. Shown is the power P as a 
function of the impurity parameter a. In all cases N = 100 
and ui — 2.5. The top panel is for C = 1.0, while the bottom 
one is for C = 0.2. The branch designation is as follows: 
a) Unstable soliton centered at the impurity (n — no), b) 
stable on-site soliton centered at n — no, c) Unstable inter- 
site soliton centered at n — no + 0.5, d) stable on-site soliton 
at n = no + 1, e) unstable inter-site soliton at n — no + 1.5, 
f) stable on-site soliton at n = no + 2, g) unstable inter- 
site soliton at n = no + 2.5, and h) stable on-site soliton at 
n = no + 3. The stable on-site mode located at the impurity, 
in the homogeneous case, disappears for a coupling value of 
C ~ 1.25 due to resonances with the phonon band. 



a = a c < 0, it also disappears together with the unsta- 
ble mode localized at n = no + 1.5 through a saddle-node 
bifurcation. If we increase again the impurity parame- 
ter, this unstable mode localized at n = no + 1.5 bifur- 
cates with the stable site mode localized at n = no + 2 
for a critical value of parameter a = a' c > through a 
saddle-node bifurcation again, and it could be possible 
to continue this bifurcation pattern until a site hq + k, 
where the value of site k increases with the value of the 
coupling C and the frequency uj parameters. This sce- 
nario is similar to the one found in previous studies with 
different kinds of impurities (40l . |4o| and appears to be 
quite general. It should be noted that when the coupling 
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parameter increases, more bifurcations take place, in a 
narrower interval of power P and impurity parameter a 
values. 

Some of the particularly interesting experimentally 
tractable suggestions that this bifurcation picture brings 
forward are the following: 

• A localized mode centered at the impurity may be 
impossible for sufficiently large attractive impuri- 
ties (because the amplitude of the mode may de- 
crease to zero), while it may be impossible to ob- 
serve also in the defocusing case due to the insta- 
bility induced by the pitchfork bifurcation with its 
neighboring inter-site configurations. 

• A localized on-site mode centered at the neighbor- 
hood of the impurity should not be possible to lo- 
calize for sufficiently large impurity strength both 
in the attractive and in the repulsive impurity case. 

We have also performed a more detailed study of the 
bifurcation between the on-site nonlinear mode centered 
at the impurity and its inter-site and one-site neighbor. 
Thus, we have determined that, for a given value of the 
coupling parameter C, the corresponding critical value of 
impurity parameter a = a c . Note that this bifurcation 
takes place only if a is negative (repulsive impurity) . In 
case of a positive (attractive impurity), the inter-site so- 
lution disappears in a saddle-node bifurcation with the 
on-site wave centered at the site next to the impurity. 
In these cases, via an analysis of invariant manifolds of 
the DNLS map, and following the method developed in 
Ref. [HI (see Appendix A) , some approximate analytical 
expressions corresponding to this bifurcation point can 
be obtained. Fig. [3] shows the comparison between the 
exact numerical and the approximate analytical results. 
In general, for a fixed value of the coupling parameter C, 
the critical value of the frequency increases with \a\. 
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FIG. 3: Bifurcation loci corresponding to the bifurcation be- 
tween the on-site localized mode at the impurity (n = no) and 
its neighbor inter-site breather (n — no + 0.5) (top panel), and 
to the bifurcation between the on-site localized mode next to 
the impurity (n = no + 1) and its neighbor inter-site breather 
(n = no + 0.5) (bottom panel), for different values of pa- 
rameter C. Dashed lines correspond to numerical results and 
continuous lines to approximate analytical calculations. 



IV. THRESHOLD FOR SOLITARY WAVE 
FORMATION 

We now examine the problem of solitary wave forma- 
tion, i.e., whether there exists a minimal, say, amplitude 
threshold for a compactum of initial data u n (0) = AS n .k 
to nucleate a localized mode. The recent work of Ref. [4§| 
suggests that a good approximation to the amplitude of a 
single-site initial condition at site k required to nucleate 
a nonlinear localized mode at that site is given by 



2 



(2C- a k )A 2 < 0. 



(13) 



In this expression, a k is the impurity parameter value at 
site k, and A the amplitude of the initial condition (48j . 

In order to study the effect of the impurity on this mag- 
nitude, we have performed numerical simulations initially 
"seeding" energy at different sites of the lattice (either at 



the impurity or at its neighbors). After a transient state, 
we have analyzed the existence of localized modes on the 
chain. To measure the localization of a state we have 
introduced the localization of an initial excitation of am- 
plitude A, L(A), as 



L(A) 



E» |yn| 2 

(£>n|) 2 ' 



(14) 



Thus, for a single excited particle we have L = 1, and 
if we have n excited particles (with the same amplitude, 
and the rest with zero amplitude), L — 1/n. In general, 
l/N<L<\. 

In Fig. [4] we summarize our numerical results and ana- 
lytical prediction. In general, when a single perturbation 
is located on the impurity, numerical and analytical re- 
sults are in good agreement (left panel of Fig. SJ. On the 
other hand, when the perturbation is located in other 
(nearby to the impurity) sites of the chain, the excita- 
tions of impurity dynamics play a significant role, and 
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FIG. 4: Localization as function of amplitude A and impurity 
parameter a for a single excitation <p n (t = 0) = A<5 ni &. The 
left panel corresponds to the excitation localized at the impu- 
rity (k = no), the center panel to the excitation localized at 
the first neighbor of the impurity (k — no + 1) , while the right 
panel to the second neighbor of the impurity (k = no + 2). 
The solid line depicts in each case the theoretical threshold 
given by Eq.[[I5]l. In all cases N = 200 and C = 1.0. 



numerical and analytical thresholds show a slight diver- 
gence (middle panel of Fig. [4j. However, when the per- 
turbation is located far enough of the impurity (that for 
the purposes of formation of a localized mode, we return 
to the limit of a "homogeneous" lattice), the effect of the 
impurity is negligible, and the threshold corresponding 
to homogeneous case is in good agreement with the nu- 
merical data, as can be appreciated in the right panel of 
Fig.rjfl 

These results also suggest an immediately testable ex- 
perimental prediction, namely that thresholds such as the 
ones reported in Ref. [49| (see also Ref. [5(| for the defo- 
cusing case) should be directly affected by the presence 
of a localized impurity. In particular, an attractive lin- 
ear impurity facilitates the formation of localized modes, 
by decreasing the threshold of their formation, while the 
opposite is true for repulsive impurities that increase the 
corresponding threshold. 



V. INTERACTION OF A MOVING LOCALIZED 
MODE WITH A SINGLE IMPURITY 

Early studies of the DNLS had shown that discrete 
solitary waves in the DNLS can propagate along the 
lattice with a relatively small loss of energy (5f|, and 
more recent work suggests that such propagating solu- 
tions might exist, at least for some range of control pa- 
rameters @, [H, nevertheless, genuinely traveling so- 
lutions are not present in the DNLS, but only in variants 
of that model (such as the ones with saturable or cubic- 
quintic nonlinearity) [58| . 



FIG. 5: Power trapping (left), reflection (center) and trans- 
mission (right) coefficients as function of impurity parameter 
a and initial thrust q. In all cases N = 1000 and C = 1. 

In this section we study the interaction of propagating 
(with only weak radiative losses) localized modes with 
the impurity. Thus, we consider a nonlinear localized 
mode, far enough from the impurity, of frequency to, and 
perturb it by adding a thrust q to a stationary breather 
, so that: 

<Pn(t = 0) = 0„e^". (15) 

This is similar in spirit to the examination of Ref. [4a |. 
although we presently examine both attractive and re- 
pulsive impurities. In the remainder of this study we 
consider ui = 2.5 and C = 1, but we have checked that a 
similar scenario emerges for other values of the frequency 

LJ. 

In general, if q is large enough, the soliton moves with a 
small loss of radiation. We have calculated, as a function 
of parameters q and a, the power and energy that remains 
trapped by the impurity, reflected and transmitted along 
the chain, and determined the corresponding coefficients 
of trapping, reflection and transmission, defined as the 
fraction of power (energy) that is trapped, reflected or 
transmitted. In Fig. [5] we summarize our results. 

We can essentially distinguish four fundamental 
regimes: 

(a) Trapping. If the parameters q and a are small 
enough, and the impurity is attractive, nearly all 
the energy remains trapped at the impurity, and 
only a small fraction of energy is lost by means of 
phonon radiation. An example of this phenomenon 
is shown in Fig. [S] (top). In this case, the cen- 
tral power (power around the impurity) before the 
collision is nearly zero. When the localized mode 
reaches the impurity, it loses power as phonon ra- 
diation and remains trapped. The analysis of the 
Fourier spectrum of this trapped breather, carried 
out after the initial decay and at an early stage of 
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FIG. 6: Trapping: Contour plot corresponding to the power 
of soliton P as function of site n and time t (top panel) and 
Fourier components of the trapped soliton calculated soon 
after the collision (bottom panel). The parameters are a = 
0.2, q — 0.3, lo — 2.5, (7=1 and the impurity is located at 
n = 0. 



the evolution, shows a frequency close to the ini- 
tial soliton frequency, as shown in Fig. [6] (bottom). 
We have observed that, in general, this frequency 
is slightly smaller than that of the incident soliton, 
and, in consequence, it has even smaller energy (in 
absolute value) and power than the corresponding 
nonlinear mode with the frequency of incident soli- 
ton. 

In this particular case, corresponding to q = 0.3 
and a = 0.2, the initial incident wave (after pertur- 
bation) has power P = 2.61 and energy E — —5.40 
and the stationary mode, trapped at the impu- 
rity, with the same frequency, has P = 2.17 and 
E = —4.73. Thus, the incident breather can acti- 
vate this nonlinear mode, and nearly all energy and 
norm remains trapped. In all simulations we have 
detected similar phenomena, as reported recently 
in a Klein- Gordon system (60. IrJlj. 

(b) Trapping and reflection. If the impurity is attrac- 
tive, but strong enough, some fraction of energy re- 



mains trapped by the impurity, but a considerable 
amount of it is reflected. The reflected excitation 
remains localized. This case is similar to the pre- 
vious one, but now the incident traveling structure 
has enough energy and norm to excite a stationary 
mode centered at the impurity, remaining localized 
and give rise to a reflected pulse. A typical case 
is shown in Fig. corresponding to q — 0.6 and 
a = 1.0. The incident wave has power and energy 
P = 2.61 and E = —4.79, and the stationary non- 
linear mode centered at the impurity, with the same 
frequency, P = 0.76 and E = —1.79. When the in- 
cident breather reaches the impurity, it excites the 
nonlinear mode, and, after losing some energy (in 
absolute value), part of it remains localized, and an- 
other part is reflected. Also, in our numerical sim- 
ulations, we have detected, as in the previous case, 
that the frequency of the remaining trapped mode 
is slightly lower than that the incident breather, so 
it has even smaller energy (in absolute value) and 
power than the corresponding nonlinear mode with 
the frequency of incident soliton. 
In general, we have found that a necessary condi- 
tion to trap energy and power by the impurity is 
the existence of a nonlinear localized mode centered 
at the impurity, with similar frequency, and energy 
(in absolute value) and power smaller than that of 
the corresponding incident soliton. 

(c) Reflection with no trapping. Here, we have to dis- 
tinguish two cases. If the impurity is repulsive, and 
q small enough, neither trapping, nor transmission 
occur. Instead, all energy is reflected, and the trav- 
eling nonlinear excitation remains localized. In this 
case, as shown in Fig. [5] (top), the incident wave 
has no energy and power to excite the localized 
mode. In a typical case, i.e., u — 2.5, q = 0.6 
and a = —0.5, the incident soliton has energy and 
power E = —4.79 and P — 2.61, and the nonlin- 
ear localized mode on the impurity with the same 
frequency E = —8.038 and P — 3.77. No trapping 
phenomenon occurs, and the pulse is reflected. 

On the other hand, if the impurity is attractive 
and strong enough, i.e. q = 0.7, u> = 2.5 and a = 
2.0, the frequency of the soliton is smaller than the 
corresponding to linear impurity mode (ul — 2.82), 
and all the energy is reflected. This is in accordance 
with the necessity of a nonlinear localized mode at 
the impurity site in order for the trapping to occur. 

(d) Transmission with no trapping. If |a| is small 
enough, and q high enough, transmission with no 
trapping occurs, as shown in Fig. [5] (bottom). 
There exists a critical value of q — q c > that, 
if q > q c , the incident soliton crosses through the 
impurity. The value of q c grows with \a\. In the 
case where q < q c , if a < 0, reflection with no 
trapping occurs, while if a > 0, trapping with no 
reflection phenomenon takes place. 
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FIG. 7: Trapping and reflection: Contour plot corresponding 
to the power of soliton P as a function of site n and time t (top 
panel) and Fourier components of the trapped soliton calcu- 
lated soon after the collision (bottom panel) . The parameters 
are a = 1.0, q = 0.6, ui — 2.5, (7 = 1 and the impurity is 
located at n — 0. 



FIG. 8: Reflection with no trapping (top panel) correspond- 
ing to parameters a = —0.5, q = 0.6 and u> — 2.5 and trans- 
mission with no trapping (bottom panel) corresponding to 
parameters a = 0.1, q — 0.7 and to = 2.5. In both cases we 
represent a contour plot corresponding to the power of soliton 
P as function of site n and time t, C = 1 and the impurity is 
located at n — 0. 



Our results related to trapping, reflection and trans- 
mission phenomena are in agreement with some results 
recently obtained, using a different approach, in a simi- 
lar system 45] . In this work, where approximate discrete 
moving solitons with fixed amplitude are generated using 
a continuous approximation, the authors study the trap- 
ping process by a linear and a nonlinear attractive impu- 
rity. In this latter framework, trapping can be explained 
by means of resonances with the linear localized mode. 
In our case, where nonlinear effects become stronger, all 
this phenomena are related with resonances with a non- 
linear localized mode. 

Finally, a very interesting phenomenon occurs when 
parameter a is repulsive and small (in absolute value) 
enough. In this case, the solitary wave can be reflected 
or transmitted depending on its velocity. Also, when it 
is reflected, our numerical tests show that its velocity 
is similar to its incident velocity. Thus, if we consider 
the soliton as a "quasiparticle" , the effect of the impu- 
rity is similar to the effect of a potential barrier. To 
determine this potential barrier for a given value of pa- 



rameter a, we have used a method similar to the one 
described by Ref. [62j • We have considered different val- 
ues of the thrust parameter q corresponding to the reflec- 
tion regime, and determine, for each value, the turning 
point, X(q). Thus the translational energy of the barrier 
for this value of q is defined as the difference between 
the energy © of the moving soliton (fT5|) and the sta- 
tionary state ([5]) of the same frequency far from the im- 
purity. It can be written as V(q) = C sm(q/2)\P(q/2)\, 
with P(q) = iY.n^ni'n+i - iK]pn-i being the lattice 
momentum, as defined in Ref. 63] . Results are shown in 
Fig. [9j which exhibits, as expected, an irregular shape, 
whose origin lies in the nonuniform behavior of the trans- 
lational velocity due to the discreteness of the system. 

On the other hand, if the parameter a is small enough, 
and positive (attractive), the solitary wave faces a poten- 
tial "well" and can be trapped if its translational energy 
is small or, if the translational energy is high enough, it 
may be transmitted, losing energy that remains trapped 
by the impurity, and decreasing its velocity. We have not 
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X 



FIG. 9: Contour plot of the phenomenon of reflection of a soli- 
ton corresponding to trust parameter q — 0.6 (top panel). Po- 
tential barrier calculated as described in text (bottom panel) . 
In both cases a — —0.2, C = 1, ui = 2.5 and the impurity is 
located at n = 0. 

found a regime with trapping and transmission as have 
also been observed in Klein-Gordon lattices [60(. 



purity significantly modifies the threshold for the forma- 
tion of localized modes, under a compactum of initial 
data. Attractive impurities favor the formation of such a 
mode under weaker excitations, while repulsive ones ne- 
cessitate an even higher amplitude threshold. Finally, we 
have examined in detail for both impurity cases (attrac- 
tive and repulsive) the interaction of the impurity with 
a moving localized mode initiated away from it. The 
principal regimes that we have identified as a function 
of the impurity strength (and sign) and initial speed are 
trapping, partial trapping and partial reflection, pure re- 
flection and pure transmission. In general, if the impu- 
rity is repulsive, and the speed small enough, the wave 
is always reflected. If the impurity strength (in absolute 
value) is small enough and the speed is high enough, then 
transmission can take place. On the other hand, if im- 
purity is attractive, trapping can occur, and if the speed 
is high enough trapping with reflection too. If impurity 
is attractive and sufficiently strong, the frequency of the 
soliton is smaller than the one corresponding to the linear 
localized impurity mode and the wave is reflected. 

There are numerous avenues that one can think of for 
further exploration of this subject. On the one hand, we 
feel that numerous among the conclusions of the present 
work including ones about the unavailability of localiza- 
tion on or at nearby sites to the impurity for sufficiently 
high strengths, or ones about the threshold for localized 
modes should be immediately experimentally testable in 
arrays of optical waveguides. On the other hand, this 
type of wave-impurity interactions have been predom- 
inantly studied in one-dimensional systems. However, 
the present availability of two-dimensional waveguide ar- 
rays renders this a very interesting system for examining 
the relevant interaction in multi-dimensional frameworks, 
even from a theoretical point of view and the examination 
of both the standing wave and of the scattering problems. 
The latter problem is currently under investigation and 
will be reported in future publications. 



VI. CONCLUSIONS 

In this work, we have revisited the long-standing theme 
of the interactions of DNLS localized modes with an im- 
purity. In particular, we have examined both the case 
of attractive and repulsive impurities and have shown 
how localized modes bifurcate out of the linear spectrum 
in the presence of the impurity. Subsequently, we have 
seen how drastically the bifurcation diagram of localized 
modes is affected by the presence of the impurity. In 
particular, we have concluded that for attractive impu- 
rities the on-site mode at the impurity eventually dis- 
appears, while for repulsive ones, it becomes unstable 
beyond a critical impurity strength. In addition, local- 
ized modes one site away from the impurity and beyond 
are also structurally affected and cannot be sustained 
under strong (either attractive or repulsive) impurities. 
Furthermore, we have seen how the presence of the im- 
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APPENDIX A: INVARIANT MANIFOLDS 
APPROXIMATION 

In this appendix we sketch the method followed in Sec- 
tion 4.1.4 of Ref. [U for determining the value of a c , i.e., 
the value of a at which the breathers centred at n = tiq 
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FIG. 10: First winding of the homoclinic tangle of the map 
(|A1|) . Dashed line corresponds to the linear transformed un- 
stable manifold when a = 0. Labels 1, 2, 3 (1', 2', 3') corre- 
sponds to fundamental solitons for a — (a/0). 



A(a) 




(A3) 



When a > 0, the unstable manifold moves downwards, 
changing the intersections between the transformed un- 
stable manifold and the stable manifold to points 1', 2' 
and 3' (see Fig. [TU]). For a — a c , both manifolds become 
tangent. Thus, for a > a c intersections 3' and 2' are 
lost, that is, for a — a c the breathers centred at n = 1 
and n = 0.5 experience a tangent bifurcation. On the 
contrary, if a < 0, intersections 1' and 2' are lost when 
\a\ > \a c \, leading to a bifurcation between the breathers 
centered at n = 0.5 and n = 0. 

A method for estimating a c (u>) is based on a simple 
approximation of W u (0). Let us consider a cubic approx- 
imation W™ pp of the local unstable manifold of Fig. [TUl 
parametrized by y = X x — c 2 x 3 , with A = A+ . The coef- 
ficient c depends on lu and C and need not be specified 
in what follows (a value of c suitable when A is large is 
computed in Ref. [HI). We have 



and n = uq + 0.5 bifurcate. 

The difference equation ([6]), for a — 0, can be recast 
as a two-dimensional real map by defining y n = <f> n and 
x n = d> n -i 0: 



V = Ao£ 



c 2 x 3 



(A4) 



on the curve A(ui, a)W a u pp , where A = A — a/C. By 
symmetry we can approximate the local stable manifold 
using the curve W£ pp parametrized by 




vD/c 



(Al) 



For uj > 2, the origin x n = y n = is hyperbolic 
and a saddle point. Consequently, there exists a one- 
dimensional stable (W s (0)) and a one-dimensional un- 
stable (W u (0)) manifolds emanating from the origin in 
two directions given by y = X±x, with 



A± = 



4C 2 



2C 



(A2) 



These manifolds intersect in general transversally, 
yielding the existence of an infinity of homoclinic orbits. 
Each of their intersections corresponds to a localized so- 
lution. Fundamental solitons (i.e. on-site and inter-site 
solitons), correspond to the primary intersections points, 
i.e. those emanating from the first homoclinic wind- 
ings. Each intersection point defines an initial condition 
(xo,yo), that is, (<p-i,(fia), and the rest of the points com- 
posing the soliton are determined by application of the 
map (|Aip and its inverse. Fig. [TU] shows an example of 
the first windings of the manifolds. Intersections corre- 
sponding to fundamental solitons are labeled as follows: 
(1) is the on-site breather centred at n = 0, (2) is the 
inter-site breather centred at n = 0.5 and (3) is the on- 
site breather centred at n = 1. 

The effect of the inhomogeneity is introduced as a lin- 
ear transformation of the unstable manifold A(a)W u (0) 
with A(a) given by: 



x = Xy 



2 „,3 



c y 



(A5) 



The curves A(a)W" p and W^ pp become tangent at (x, y) 
when in addition 



(A-3c 2 x 2 )(A -3c 2 y 2 ) = 1. 



(A6) 



In order to compute a c as a function of u), or, equiv- 
alently, the corresponding value of Ao as a function of 
A, one has to solve the nonlinear system (|A4[) ~ (|A6[) with 
respect to x, y, Ao, which yields a solution depending on 
A. Instead of using A it is practical to parametrize the 
solutions by t = y/x. This yields 



1 



iV2 



t 3 ' 



,V2 



t 



y 



3 

Ao = ^ 



1 

2^ : 



(*• 



1 



t 3 ' 



A = h- - t 



2t 2 

Since A + A -1 = oj/C it follows that 
t 4 - 2Xt + 3 = 0, 



2 v t ' 



(A7) 



(A8) 



Given a value of to, one can approximate a c by the value 
of a given by equations (|A7[) - (|A8[) . In particular, (|A7|) 
has two real positive solutions (one larger than 1, and 
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another smaller than 1), and two complex conjugated 
solutions. The solution with t > 1 (t < 1) leads to a c > 
(a c < 0) and, subsequently, approximates the tangent 
bifurcation values when the breathers at n = 0.5 and 
n = 1 (n = 0) collides. 



Despite it gives precise numerical results in a certain 
parameter range, the approximation (|A7[) - (|A8[) is not al- 
ways valid. Indeed, the parameter regime ui < 5C/2 is 
not described within this approximation [651 ]. 
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